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Abstract 

The  effects  of  a  uniform  magnetic  field  normal  to  the  shock  front  upon  shock  states 
and  expansion  wave  properties  are  theoretically  and  experimentally  investigated  for 
large  magnetic -drive  fields.  One -dimensional  theory  and  several  simplifying  assump¬ 
tions  are  employed  to  predict  the  important  effects,  and  a  simple  method  of  visualizing 
the  possible  shock  states  is  obtained.  The  restrictions  upon  shock-wave  formatim 
imposed  by  various  types  of  instabilities  are  considered.  A  series  of  experiments  has 
been  performed  in  a  magnetically  driven  shock  tube  of  annular  cross  section  under  c(»i- 
ditions  approaching  the  assumptions.  Some  difficulty  was  encountered  in  producing  well- 
formed  shocks,  and  a  tentative  link  with  the  drive-current  breakdown  characteristics  is 
established.  Observations  of  well-formed  shocks  indicate  close  agreement  with  shock- 
transition  theory,  even  when  the  axial  magnetic  field  exceeds  the  switch-on  value, 
although  the  expansion  wave  properties  exhibit  large  departures  from  the  theory.  Some 
estimates  of  flow  adjustments  caused  by  radius  effects  and  leakage  through  the  boundary 
layer  have  been  made,  and  are  found  to  contribute  to  an  understanding  of  the  experi¬ 
mental  results. 
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I.  INTRODUCTION 


The  magnetically  driven  shock  tube  is  cap>able  of  producing  extremely  large  veloc¬ 
ities  and  high  temperatures,  as  well  as  complicated  magnetohydrodynamic  shock  waves, 
and  consequently  is  the  present  concern  of  intensive  theoretical  and  experimental  inves¬ 
tigation.*"^  Because  magnetohydrodynamic  shock  waves  couple  the  laws  of  gasdynamics 
with  those  of  electromagnetism,  the  geometry  of  the  shock  tube  must  be  carefully  chosen 
to  make  a  theoretical  analysis  possible.  The  magnetic  annular  shock  tube  (MAST)  devel¬ 
oped  1^  the  Avco-Everett  Research  Laboratory  produces  shock  waves  that  approach  the 

1  -3 

one -dimensional  model,  for  which  case  a  complete  analysis  is  available.  Such  a 
device  was  constructed  in  the  Magnetogasdynamics  Laboratory  of  the  Research  Labora¬ 
tory  of  Electronics,  M.  I.  T. ,  and  an  experimental  program  was  initiated  to  gain  an 
understanding  of  its  actual  performance  and  to  explain  any  major  deviations  from  the 
theory. 

Long  before  the  shock  tube  was  functioning,  the  effect  of  a  magnetic  field  normal  to 

the  shock  front  (axial  magnetic  field)  was  chosen  for  investigation.  Some  of  the  most 

interesting  magnetohydrodynamic  shock  phenomena  can  occur  in  the  presence  of  a 

normal  magnetic  field.  These  phenomena  include  increased  shock-front  velocities, 

spreading  of  the  expansion  wave,  and  switch-on  shock-wave  transitions  in  which  part 

of  the  drive  current  flows  inside  the  shock  front.  *  The  question  has  also  been  raised 

recently  about  what  type  of  shock  transition,  if  any,  is  stable  when  the  axial  magnetic 

2 

field  exceeds  the  switch-on  value. 

Since  future  shock-tube  experiments  are  likely  to  incorporate  axial  magnetic  fields 
for  such  purposes  as  acceleration  of  the  shock,  stabilization  of  the  boundary  layer,  or 
containment  of  the  plasma,  it  is  necessary  to  determine  whether  these  fields  have  any 
unexpected  or  undesirable  effects  upon  magnetically  driven  shock-tube  performance. 

The  analyses  and  experiments  reported  here  were  conducted  to  investigate  the  effects 
and  answer  the  questions  discussed  above. 

Section  II  of  this  report  is  devoted  to  the  theory  of  magnetohydrodynamic  shock  waves 
and  expansion  waves  and  includes  all  information  necessary  for  predicting  shock-tube 
performance  based  upon  one -dimensional  theory.  Section  Ill  contains  a  detailed  descrip¬ 
tion  of  the  experimental  apparatus,  as  well  as  of  the  experimental  procedure.  Section  IV 
presents  the  experimental  data  in  both  raw  and  reduced  forms,  relates  it  to  the  theory, 
and  mentions  some  of  the  important  observations  on  the  functioning  of  the  experimental 
equipment.  Section  V  completes  the  main  part  of  the  report  by  giving  conclusions  and 
recommendations . 

Appendix  I  is  devoted  to  some  estimates  of  deviations  from  one-dimensic»ial 
theory  which  arise  from  radius  effects  and  boundary -layer  leakage.  Appendix  II  con¬ 
tains  absolute  light  intensity  information  for  use  in  measuring  emitted  exptirimental 
light  and  predicting  its  magnitude.  Appendix  III  gives  a  complete  list  of  the  auxiliary 
equipment  required  for  this  experiment. 
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The  naks  system  is  employed  generally  throughout  the  report,  except  for  some  com¬ 
puted  values  that  are  stated  in  more  convenient  units. 
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II.  ONE -DIMENSIONAL  MAGNETOHYDRODYNAMIC  SHOCK-WAVE  THEORY 


Simple  magnetohydrodynamic  shock-wave  theory  will  now  be  outlined,  and  equations 
derived  which  are  required  for  comparison  with  the  experimental  results.  Some  empha¬ 
sis  is  placed  upon  making  the  phenomena  intuitively  accessible.  In  order  to  emphasize 
the  relationship  between  the  theoretical  model  and  the  experimental  situation,  a  brief 
description  of  the  experiments  will  be  given. 

The  experiments  were  designed  to  furnish  shock-velocity  and  compression-ratio  data 
for  high-speed  shock  waves  (M  »1)  under  the  influence  of  constant  magnetic  fields  perpen¬ 
dicular  to  the  shock  front  (axial  fields).  The  magnitudes  of  the  axial  magnetic  field  (B^) 
and  experimental  gas  density  (Pj)  were  chosen  before  each  experiment  from  a  wide  range 
of  possible  values.  The  shock-tube  geometry  was  such  that  one-dimensional  theory  is  a 
reasonable  first-order  approximation,  all  properties  varying  only  in  the  axial  (x)  direction. 

The  magnetically  driven  shock-wave  phenomenon  consists  of  two  physically  dependent 
but  conceptually  separable  phenomena:  the  shock  layer  and  the  expansion  wave  (containing 
the  drive  current).  Figure  1  qualitatively  indicates  the  impx)rtant  features  of  the  one- 
dimensional  shock  wave,  as  well  as  the  coordinate  system  used  throughout  this  repjort. 
The  flow  is  from  -x  to  +x(u>0).  Even  though  the  annular  shock  tube  has  cylindrical  sym¬ 
metry,  Cartesian  coordinate  system  nomenclature  is  used  as  a  reminder  of  the  assump¬ 
tions  made  in  the  analysis.  Careful  distinction  is  maintained  between  the  expaerimental 
gas  (that  is,  the  gas  present  in  the  shock  tube  before  the  experiment  commences  amd 
upstream  of  the  shock  front  after  the  expjeriment  commences)  and  the  test  gas  (that  is, 
the  gas  present  between  the  shock  front  and  the  upstream  edge  of  the  expansion  wave). 
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Fig.  1.  Coordinate  system  for  one -dimensional  shock  theory. 


2.  1  THE  SHOCK  LAYER 
a.  Stationary  States 


The  fluid  properties  upstream  and  downstream  of  the  shock  layer  are  stationary  and 
are  related  through  the  conservation  equations.  These  propierties  are  independent  of 
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all  transport  coefficients,  which  only  control  the  thickness  of  the  shock  layer, 
has  reduced  the  conservation  equations  to  the  form  required  for  this  ainalysis. 

Shercliff 

Assuming 

(a)  steady  flow  (shock-stationary  reference  frame). 

(b)  variation  in  the  x -direction  only  (one 

(c)  isotropic  gas  pressure, 

(d)  a  perfect  gas  as  the  working  fluid, 

(e)  v  #  0, 

(f) 

-dimensional  shock). 

Shercliff  finds 

F  B 

y  X 

B  (u-u  )  =  — Ti —  =  constant 
y  0  o 

(transverse  momentum) 

(1) 

B^ 

y 

p  +  Gu  +  =  F^  =  constaint 

(Rayleigh) 

(2) 

^  +  T  Y  +  “  constant 

(Fanno) 

(3) 

B^  a  constant 

(4) 

G  a  pu  a  constant 

(5) 

u  a  =  constant 

o  pG 

(6) 

B  a  W  a  0 
z 

(7) 

Ja  E  a  V  XB  a  0. 

(8) 

These  relations  are  also  valid  within  the  shock  layer,  except  that  j^  =  <r(uBy-vB^)  #  0. 
When  B  =0  the  relations  form  the  usual  gasdynamic  shock  relations,  even  if  B  is 

y  B  * 

allowed  to  approach  zero  When  u  =  u  ,  Eq.  6  shows  that  u  =  — —  =  b  . 

In  order  to  solve  these  equations  outside  the  shock  layer  and  obtain  the  possible 
stationary  states  for  the  fluid,  it  is  necessary  to  specify  the  values  of  seven  of  the  func¬ 
tions.  Since  B  ,  p,  k,  and  p  are  known  in  advance  of  an  experiment,  the  knowledge  of 
any  additional  three  (for  example,  F  ,  F  ,  and  H,  or  u,  p,  and  B  )  is  required  for  a 

2  *  y  y 

solution.  Following  Anderson,  a  knowledge  of  F  ,  F  ,  and  H  is  assumed,  and  the 

X  y 

solution  of  Eqs.  1-3  is  represented  in  three  dimensions.  Greatest  clarification  of  the 

2 

experimental  problem  follows  from  the  choice  of  u-p-B  space  for  this  representation 
(Fig.  2). 
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(0) 


(b) 


Fig.  2.  Graphical  representation  of  stationary  shock  states.  (a)B  ,^0.  {b)B 


« 


(i)  Transverse  momentum  surface.  This  is  a  cylinder,  with  generators  parallel 

2 

to  the  p-axis.  In  the  u-B  plane  it  resembles  a  hyperbola,  the  asymptotes  being  u  b  u 

2  y  2  ° 

and  By  a  o.  When  Fy  b  0  it  degenerates  to  the  planes  u  »  and  By  b  o.  Since  Eq.  1 
has  been  squared  to  suit  this  representation,  it  must  be  remembered  that  By  changes 
sign  as  u  passes  through  the  value  u^. 

(ii)  Rayleigh  surface.  This  is  a  plane  intersecting  the  coordinate  axes  at  u  = 
p  =  Fjj,  and  By  =  ZpFj^. 

(iii)  Fanno  surface.  This  is  the  most  complex  of  the  three,  and  it  will  be  noted  only 

that  this  surface  is  asymptotic  to  u  s  0,  it  intersects  the  u-axis  at  its  maximum  u-value 
of  and  it  has  negative  slope  and  positive  curvature  in  any  plane  of  constant  By. 

Referring  to  Fig.  2a,  it  is  evident  that  there  are  no  more  than  four  stationary  states 

for  a  given  set  of  B  ,  p,  k,  p,  F  ,  F  ,  and  H.  They  are  numbered  in  the  conventional 
X  ^  X  y 

manner,  which  may  be  shown  to  be  the  direction  of  decreasing  kinetic  energy  and  veloc¬ 
ity,  and  increasing  density,  pressure,  and  entropy.  The  second  law  of  thermodynamics 
restricts  the  allowable  combinations  of  the  four  equilibrium  states  to  six,  namely,  1-2, 
1-3,  1-4,  2-3,  2-4,  and  3-4. 

Decreasing  Fy  from  the  value  in  Fig.  2a  causes  state  2  to  approach  state  3  until, 

an 


If 


when  states  2  and  3  coincide,  F  =  0  and  u,  =  u,  =  b  (Fig.  2b).  This  is  called 

y  3  X 

Alfv4n  shock  and  involves  only  a  change  in  sign  of  B  .  It  follows  that  u,  >  b  ^  u-. 

y  fa  X  j 

Fy  is  Increased  over  the  value  of  Fig.  2a  until  either  states  1  and  2  or  states  3  and  4 
coincide,  Eqs.  1-3  reveal  that 


u^  -  u^(a^+b^)  +  b^a^  =  o  (magnetoacoustic  wave-speed  equation). 


(9) 


Here,  either  Uj  =  u^  =  or  u^  s  a  Consequently,  Uj  >  Cj^  >  u^.  and  >  c^^ 

^  u^.  The  six  allowable  shock  transitions  of  the  diagram  may  now  be  characterized; 

1-2  (fast  shock): 

3-4  (slow  shock); 


increases,  u^  ^  ?  u^  ^ 


iByl 

B  1  decreases,  b  ^u,  Sc  ^u.; 

I  y'  '  X  3  sx  4 

1-3,  1-4,  2-3,  2-4  (intermediate  shocks);  By  and  (u-u^)  change  sign. 


In  this  experiment  the  observed  shock  velocity  (u^)  exceeds  Cj^  in  the  experimental 

gas,  so  that  u^  corresponds  to  Uj,  and  the  state  of  the  experimental  gas  corresponds 

to  state  1  of  the  theory.  Since  B  ,  b  0  and  B  ^0,  F  =  0,  and  the  transverse  momen- 

y  1  X  y 

turn  surface  takes  on  the  degenerate  form  discussed  above  and  shown  in  Fig.  2b.  Allowing 
only  B^  (that  is,  u^)  to  vary  in  Fig.  Zb,  we  see  that  it  is  obvious  that  states  2  and  3  do 
not  exist  when  u^  «  u^  or  S  Uj,  so  that  the  only  possible  shock  trsmsition  is  the  gas- 
dynamic  shock  (1-4).  However,  when  u^  «  u^  «  Uj  there  exist  four  equilibrium  states, 
two  of  which  have  By  0.  Because  of  the  unusual  behavior  of  By  in  this  case,  the  1-2 
and  1-3  shocks  are  called  "switch-on"  shocks,  and  the  2-4  and  3-4  shocks  are  called 
•switch-off"  shocks. 

2 

The  algebra  may  be  simplified  by  observing  that  Mj  »  1,  whence  PjUj  »  pj,  and  Pj 
may  be  neglected.  It  is  therefore  possible  to  specify  the  three  required  quantities  for 
any  given  experiment  ^F^  s  p  ^Uj  b  PjU^,  Fy  =  0,  and  HBU^/Z^and  solve  the  conservation 
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equations  for  possible  states  downstream  of  the  shock  layer.  The  following  solutions 
are  easily  obtained: 

Fast  or  gasdynamlc  shock 


y4 


v,=  0 


P4 


(10) 

(11) 


(k^+1)^  ^4 


(12) 

(13) 


Intermediate  or  swltch-on  shock 


The  swltch-on  shock  can  only  occur  when 


O  4 


u 

s 


(19) 


or 


(20) 


Since  the  internal  energy  upstream  of  the  shock  Is  negligible,  these  solutions  are 

only  functions  of  the  gas  properties  downstream  of  the  shock.  The  stationary-state 

numbering  system  described  above  will  be  used  throughout  this  report. 

2 

The  u-p-By  diagram  may  also  be  used  to  represent  the  effect  of  simple  momentum 
(AF),  or  of  energy  (AH)  sources  or  sinks  in  cases  in  which  an  algebraic  solution  would 
be  overly  complicated.  For  example,  if  we  assume  that  Pj  =  =  0  and  p^  and  Pj 
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are  known,  the  effect  of  a  specific  enthalpy  source  of  strength  AH  may  be  written 


where  Uj  and  are  the  perturbed  values, 
b.  Stability 

Finding  the  possible  equilibrium  states  for  a  given  flow  is  a  simple  matter  in  com¬ 
parison  with  determining  shock-layer  stability.  In  an  extensive  treatment  of  the  subject, 
2 

Anderson  examines  under  similar  assumptions  the  six  allowable  shock  transitions  and 
draws  the  following  general  conclusions: 

(i)  When  u_  >  u,  or  u  <  u.  the  1-4  shock  is  stable. 

(ii)  When  ^  >  u^  only  the  1-2  and  3-4  shock  layers  are  stable  in  the  evolution¬ 

ary  sense.  That  is,  only  these  two  shock  layers  emit  precisely  the  number  of  waves 
required  to  conserve  all  necessary  physical  quantities  at  the  shock  layer.  Hence,  they 
can  evolve  from  the  accumulation  of  small  waves  and  can  adjust  to  unsteady  flows. 

(iii)  When  u^b^  (for  example,  behind  a  switch-on  shock),  the  shock  layer  is  unstable 
in  the  sense  that  the  equilibrium  amplitude  of  the  emerging  waves  is  unbounded  when  a 
monochromatic  wave  train  of  any  finite  amplitude  overtakes  the  layer. 

For  B  0  and  small  values  of  the  axial  magnetic  field,  u^  <  u^,  and  the  stability 
arguments  (i)  agree  with  the  conclusion  that  the  gasdynamic  shock  (1-4)  must  appear. 
For  Byj  B  0  and  large  values  of  the  axial  magnetic  field,  Uj  >  u^  >  u^,  and  the  switch-on 
shock  (1-4)  should  appear  (stability  argument  (ii)).  However,  the  switch -on  shock  may 
disintegrate  (stability  argument  (iii))  under  the  influence  of  the  small  waves  that  are 
unavoidable  in  any  real  experiment.  Because  the  entropy  increase  is  maximum  for  the 
gasdynamic  shock,  the  possibility  of  observing  this  transition  cannot  be  discounted. 
Consequently,  there  is  room  for  speculation  as  to  what  type  of  shock  transition  will  be 
observed  when  Uj  >  u^  >  u^,  if  indeed  any  shock  wave  is  observed. 

2.  2  THE  EXPANSION  WAVE 

The  expansion  wave  may  be  regarded  as  a  complex  piston  that  drives  the  shock  wave 
into  the  experimental  gas.  The  drive  pressure  is  due  to  the  interaction  of  the  distrib¬ 
uted  drive  current  (j^)  with  its  own  azimuthal  magnetic  field  (B  ),  resulting  in  a  volume 
force  of  magnitude  j^B^  in  the  negative  x-direction.  The  net  result  is  called  magnetic 
driving.  Since  the  drive  current  may  move  through  the  compressed  gas  both  by  wave 
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propagation  and  diffusion,  the  resulting  mathematical  model  is  intractable.  Fortunately, 
shock-tube  magnetic  Reynolds  numbers  are  considerably  larger  than  unity  (~10  for  this 
experiment);  thus  the  effects  of  diffusion  can  be  neglected.  This  is  entirely  equivalent 
to  the  assumption  of  infinite  electrical  conductivity. 

Under  the  following  assumpticms; 

(a)  variation  in  the  x-direction  only, 

(b)  isotropic  gas  pressure, 

(c)  a  perfect  gas  as  the  working  fluid,  k  s  J-, 

(d)  and 

(e)  a  step  function  drive  current  confined  to  a  sheet  at  t  s  0,  Kemp  and  Petschek 
computed  all  necessary  expansion-wave  properties.  A  description  of  the  important 
results  for  the  s  o  case  follows. 

Any  motion  of  the  drive  current  into  the  test  gas  must  now  depend  upon  wave  propa¬ 
gation.  When  #  0  two  longitudinal  magnetoacoustic  waves  arise,  their  velocities  rel¬ 
ative  to  the  fluid  being  and  ±®g(3iow)'  inspection  of  the  property  changes 

2 

across  these  waves  indicates  that  p  and  B^  change  in  the  same  direction  across  the  fast 
wave,  and  in  opposite  directions  across  the  slow  wave.  Since  the  expansion  wave  sep- 

b2 

y 

arates  two  regions  for  which  p  +^p^is  roughly  conserved,  only  slow  waves  may  be  pres¬ 


ent  in  the  expansion  wave. 

When  Uq  <  u^,  B^^  =  0,  e  a  >  u^,  and  c^^^  b  b^^.  When  *  0,  c^^  =  0,  and  the 
expansion  wave  cannot  be  spread  in  the  x-direction  by  the  slow  wave.  There  is  no 
physical  reason  for  spreading  under  this  condition,  however,  for  every  element  of  the 
expansion -wave  volume  force  is  directed  toward  the  shock  front,  tending  simply  to  main¬ 
tain  the  shock-front  velocity.  The  drive  current  remains  confined  to  a  single  sheet. 

When  Bj^  #  0  an  azimuthal  expansion-wave  volume  force  of  magnitude  arises. 

The  azimuthal  force  is  balanced  (»ily  by  an  inertial  force,  resulting  in  an  azimuthal 
acceleration  of  the  fluid.  Since  the  field  lines  behave  as  though  they  were  attached  to 
the  fluid  particles  (assumption  (d)),  the  azimuthal  motion  of  the  particles  within  the 
exp)ansion  wave  causes  a  distortion  of  the  field  lines.  This  distortion  carries  the  drive 
current  upstream  and  downstream,  spreading  the  drive  current  over  a  greater  quantity 
of  fluid  and  thereby  reducing  the  azimuthal  accelerations.  The  front  edge  of  tli?  expan¬ 
sion  wave  now  moves  relative  to  the  fluid  with  the  speed  *  b^^  <  u^,  and  the  fluid 
that  has  entered  the  expansion  wave  accelerates  in  the  -f x-direction.  The  presence  of 
the  axial  magnetic  field  has  changed  the  basic  character  of  the  piston.  Since  Newton's 
laws  must  be  valid  for  this  system  and  momentum  is  being  lost  through  the  spreading 
of  the  expansion  wave,  the  shock  must  accelerate  (p^  increases).  The  computations  indi¬ 
cate  that  this  acceleration  is  slight  for  the  range  of  values  of  axial  magnetic  field  avail¬ 
able  in  this  experiment  (Fig.  20).  The  computations  also  indicate  that  the  drive  current 
is  largely  concentrated  at  the  front  of  the  expansion  wave  (that  is,  just  downstream  of 


9 


the  test  gas),  which  means  that  the  Huid  properties  change  rapidly  there.  The  decrease 
of  density  and  temperature  through  the  expansion  wave  makes  assumptions  (b)  and  (d) 
less  reasonable  as  x  increases.  The  transverse  velocity  (v)  at  the  downstream  edge 
of  the  exp>an8ion  wave  is  found  to  be  of  the  order  of  u  j . 

As  is  increased  with  u^  <  u^,  b^^  increases  proportionately  (Eq.  11),  and  the 
front  edge  of  the  expansion  wave  consumes  more  of  the  test  gas.  When  becomes 
large  enough  to  make  u^  =  u^,  the  shock  front  is  no  longer  separate  from  the  expansion 
wave.  This  is  seen  to  be  the  switch-on  condition  (Eq.  19). 

For  even  greater  values  of  B^,  shock-transition  algebra  indicates  several  possible 
alternatives  for  the  test-gas  equilibrium  condition.  Expansion  wave  theory  shows  that 


the  gasdynamic  shock  (1-4)  will  have  no  test  gas  because  u^  <  b^^  =  If  the  switch- 

on  shock  appears,  u^  =  b^2  (EQ-  b^^  ^  *^3x2  expansion  wave  sep¬ 


arates  from  the  shock  layer  and  recedes  downstream.  Currents  continue  to  flow  in  the 

2 

shock  front,  carried  there  by  the  fast  wave  (both  p  and  B  increase).  As  B  is  increased, 

.  \  y  ^ 

c  p/u-  decreases  continuously,  and  the  amount  of  test  gas  increases. 


2.3  CYLINDRICAL  GEOMETRY 


The  one -dimensional  theory  is  valid  for  a  magnetically  driven  shock  moving  axially 
between  two  concentric  cylinders  of  radii  r^  and  r.,  provided  that 

(a)  the  drive  current  has  azimuthal  uniformity, 

(b)  the  drive  current  source  and  sink  are  downstream  of  the  entire  expansion  wave. 


and 

(c) 


2(r  -r.) 

-Q _ 11. 

(Vri) 


0. 


Only  (c)  presents  an  experimental  problem  because  r^  is  limited  by  the  size  of  the 
equipment,  and  r^  -  must  be  large  enough  to  provide  sufficient  test  gas  for  exp>eri- 

2(r  -r.) 

mental  observation.  The  value  of  .  >  was  reduced  to  approximately  l/3  for  this 

'  o  i’  . 

experiment,  which  value,  according  to  Patrick,  is  small  enough  to  make  one -dimensional 
theory  approximately  correct. 


2.4  QUASI-ST ATIONARY  ASSUMPTION 


The  one -dimensional  theory  is  strictly  valid  only  for  the  case  of  constant  total  drive 
current  (J).  Because  the  experiment  fails  to  provide  a  constant  J,  the  establishment 
of  equilibrium  must  require  a  finite  amount  of  time.  All  shock  properties  will  be  com¬ 
puted  from  the  one-dimensional  theory  by  using  the  instantaneous  value  of  J.  In  these 
experiments  the  current  is  fairly  constant  for  approximately  3  psec,  and  the  time 
required  for  the  acoustic  wave  to  travel  from  the  expansion  wave  to  the  shock  front  is 
approximately  0.  5  psec  (fully  ionized  hydrogen  at  3  X  10^  K).  Consequently,  the  quasi¬ 
stationary  assumption  seems  justified. 
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III.  DESCRIPTION  OF  THE  EXPERIMENT 


3.  1  EQUIPMENT 

Tiie  shock-tube  consists  essentially  of  two  concentric  stainless-steel  tubes  that 
have  been  blanked  off  at  both  ends  to  allow  the  removal  of  air  from  the  annular  volume 
separating  them  (Fig.  3).  At  one  end  there  Is  a  connection  to  the  pumping  system.  At 
the  opposite  end  are  two  tungsten-silver  (Elkonite)  electrode  rings  (Figs.  3  and  4)  that 
carry  the  drive  current  to  and  from  the  experimental  gas.  They  provide  the  shortest 
radial  path  between  the  two  cylinders,  and  thereby  induce  the  initial  breakdown  of  the 
experimental  gas  to  take  place  at  a  known  axial  position.  Since  the  shock  and  expansion 
wave  move  away  from  the  electrode  end  and  toward  the  pump  end,  the  pump  end  is 
referred  to  as  "upstream,"  and  the  electrode  ring  end  as  "downstream,"  the  use  of 
these  terms  being  consistent  with  the  nomenclature  of  the  theory. 

Slightly  downstream  of  the  electrode  ring  the  experimental  gas  was  introduced 
through  a  controlled-leak  valve  (Fig.  4)  and  flowed  slowly  (~  1  m/s)  toward  the  pump 
end.  Hydrogen  was  chosen  as  the  experimental  gas  for  two  important  reasons.  First, 
the  ionization  potential  for  the  atomic  hydrogen  electron  is  less  than  13.  6  ev. 
Consequently,  for  the  anticipated  shock  velocities  (~  10  cm/psec),  hydrogen  is  fully 
dissociated  and  ionized,^  and  no  complicated  secondary  ionizations  can  arise.  Second, 


Fig.  3.  Schematic  representation  of  the  shock  tube  and  auxiliary  equipment. 
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Fig.  4.  Cross-section  view  of  the  ring-electrode  end  of  the  shock  tube. 


the  low  density  of  hydrogen  makes  large  shock-front  velocities  possible,  with  the  result 
that  the  displacement  of  the  shock  front  during  the  short  period  of  the  experiment 
(~5  psec)is  comparable  to  the  length  of  the  apparatus.  Note  that  the  period  of  the  drive 
current  is  largely  a  function  of  the  inductance  and  capacitance  of  the  external  circuit 
because  the  inductance  does  not  increase  substantially  as  the  drive  current  moves 
upstream. 

The  range  of  experimental  gas  pressures  (pj)  was  80  p  to  200  p  (±  5  per  cent),  and 
measurements  were  made  at  a  point  near  the  pump  inlet  (Fig.  3).  Viscous-flow  calcu¬ 
lations  indicated  a  pressure  drop  of  less  than  1  p  from  the  electrode  ring  to  the  pump 
inlet  for  these  conditions,  with  the  result  that  axial  variation  of  experimental  gas 
properties  could  be  neglected.  Calculations  based  upon  the  minimum  attainable 
pressure  (with  the  experimental  gas  valve  closed)  and  available  pump  characteristics 
indicated  that  impurities  in  the  experimental  gas  totaled  less  than  1  per  cent  by  volume. 

The  electrical  currents  required  for  each  experiment  were  stored  in  capacitor 
banks  and  triggered  in  sequence  by  high-voltage  pulses  applied  to  open- spark  gaps 
(Fig.  5). 

The  drive-current  capacitor  bank  contained  172  pf  (Fig.  7)  and  supplied  an  approxi¬ 
mately  sinusoidal  current,  with  a  measured  quarter-cycle  time  of  4.75  psec  and  a 
measured  peak  value  of  49  amps/volt.  Since  the  initial  drive-current  capacitor  bank 
voltage  for  all  experiments  was  4. 5  kv,  the  peak  drive  current  was  therefore 

5 

2.2  X  10  amps  (±5  per  cent).  The  inner  electrode  was  positive  with  respect  to  the  outer 
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Fig.  5.  Electronic  sequence  and  circuit  diagram. 


electrode.  Although  no  direct  experimental  evidence  has  resulted  from  this  experi¬ 
ment,  it  is  possible  that  the  drive-current  breakdown  behavior  is  strongly  influenced 
by  the  mechanical  details  of  the  downstream  end  of  the  shock  tube,  particularly  the 
geometry  and  composition  of  the  drive-current  path  (Fig.  4).  A  lucite  endplece  was 
added  to  the  outer  cylinder  at  the  extreme  upstream  end  of  the  shock  tube  (Fig.  7)  to 
prevent  the  drive  current  from  shorting  out  directly  across  the  inner  face  of  the  lucite 
end  plate. 

The  breakdown  of  the  experimental  gas  at  the  electrode  ring  was  aided  by  the  free 
electrons  resulting  from  a  small  preionization  current  triggered  ~3  psec  in  advance 
of  the  drive  current  (Figs.  5  and  8).  The  preionization  current  had  a  peak  value  of 
approximately  10  amps  and  a  frequency  of  approximately  1  me. 
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sensing  head  (above);  stainless-steel  outer  cylinder;  lucite  axial  field  solenoid  support;  loop  (center); 

*  loop  (slightly  below  shock-tube  axis);  optical  openings  to  test  section;  preionization  current  capacitor 
power  supply  (above);  drive-current  capacitor  bank,  with  spark  gap  and  current  plates  (above);  hydrogen 
bottle  (above);  and  axial-field  current  capacitor  bank  and  spark  gap  (below). 
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Fig.  8.  Detail  facing  photomultiplier  side.  Left  to  right,  below:  preionization  current  circuit  and 
spark  gap;  mechanical  forepump;  diffusion  pump;  molecular  backstreaming  trap;  and 
diffusion  pump  valve.  Left  to  right,  above:  two  photomultiplier  units;  vacuum  sensing 
head;  standard  incandescent  lamp;  two  photomultiplier  power  supplies;  and  vacuum  gauge. 


The  axial  magnetic  field  was  produced  by  means  of  a  solenoid  that  was  external  to 
the  outer  cylinder  (Figs.  3,  4,  5,  and  7).  The  solenoid  was  also  counterwound  inside  the 
inner  cylinder  in  order  to  reduce  the  inductance  of  the  circuit.  The  axial-field  current 
capacitor  bank  (Fig.  7)  had  48  pf.  and  the  current  produced  an  observed  mean  peak  axial 
field  of  0.  034  w/m  kv  (±10  per  cent)  with  a  quarter-cycle  time  of  220  psec.  (This  was 
observed  to  be  longer  than  the  time  required  for  the  axial  field  to  penetrate  the  steel 
cylinders.)  The  axial-field  current  was  triggered  approximately  215  psec  in  advance  of 
the  drive  current,  and  the  resulting  variation  of  the  axial  field  during  the  period  of  the 
drive  current  was  observed  to  be  less  than  2  f)er  cent  of  the  mean  peak  axial-field 
value.  The  eixial  field  was  carefully  probed  for  axial  and  radial  variation  of  magnitude 
and  found  to  be  within  10  per  cent  of  the  mean  peak  field  value  for  all  points  more  than 
15  cm  upstream  of  the  electrode  ring.  Fringing  of  the  axial  field  resulted  in  larger 
variations  downstream  of  this  point. 

The  preionization  radial  electric  field  and  the  axial  magnetic  field  act  together  to 
make  the  free  electrons  of  the  preionization  current  move  in  the  azimuthal  direction 
(Hall  current), and  thereby  smooth  out  any  nonuniformities  in  the  azimuthal  distri¬ 
bution  of  the  preionization  current.  This  effect  varies  inversely  with  the  magnitude  of 
the  axial  magnetic  field  and  encourages  a  uniform  initial  azimuthal  distribution  of  the 
drive  current.  In  the  absence  of  an  axial  magnetic  field,  there  is  no  azimuthal  particle 
motion. 

Small  quartz  portholes,  located  along  both  sides  of  the  outer  cylinder  in  the  hori¬ 
zontal  plane  containing  the  shock-tube  axis,  were  provided  to  allow  visual  observation 
of  the  test  gas  (Figs.  3  and  7).  The  portholes  were  located  at  15.  3-cm  intervals  along 
the  axis,  the  first  on  each  side  being  9.  0  cm  upstream  of  the  electrode  ring.  The  inner 
surfaces  of  the  portholes  were  designed  to  be  flush  with  the  inner  surface  of  the  outer 
cylinder  in  order  to  prevent  them  from  interfering  with  the  test-gas  flow  pattern.  The 
inner  cylinder  was  left  intact.  The  outer  axial  magnetic  field  solenoid  was  wound  so 
that  the  light  emitted  by  the  hot  test  gas  could  radiate  freely  within  a  cone  centered 
about  the  radius  through  the  portholes. 

The  intensity  of  the  emitted  light  was  continuously  registered  by  two  sensitive, 
rapid-response  photomultipliers  that  were  focused  through  quartz  portholes  on  a  radial 
piortion  of  the  annular  experimental  spwce  (Figs.  3  and  8). 

Well-formed  shocks  are  recognized  by  their  discontinuous  nature.  The  ideal  tran¬ 
sition  time  in  a  coordinate  system  fixed  to  the  shock  tube  is  the  shock-layer  thickness 
divided  by  the  shock  velocity. 

The  expected  ideal  transition  time  for  this  experiment  was  ~  0.  05  psec,  but  the 
observed  value  was  much  greater,  for  the  following  reasons: 

1.  Even  if  a  step-function  light  signal  were  applied  to  a  photomultiplier,  the 
output  signal  would  be  delayed  by  the  photomultiplier  output  circuit  rise  time.  The 
rise  time  for  the  photomultiplier  circuits  used  in  this  experiment  was  observed  to  be 
~  0.  1  psec. 
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2.  In  order  to  obtain  a  useful  signal  from  the  photomultipliers  during  the  experi¬ 
ments,  it  was  necessary  to  collect  light  from  a  finite  volume  of  test  gas.  The  time 
required  for  the  shock  front  to  traverse  this  volume  further  increased  the  total  rise 
time.  The  traversal  time  was  reduced  to  ~  0.  01  psec  by  placing  small  apertures  over 
the  quartz  portholes. 

3.  The  theory  presented  in  Appendix  I  indicates  that  the  shock  front  was  not  located 
in  a  plane  normal  to  the  shock-tube  axis,  but  was  tilted  downstream  with  increasing 
radius.  This  tilt  could  have  added  as  much  as  0.  1  psec  to  the  total  rise  time,  but  the 
effect  was  reduced  by  tilting  the  axis  of  the  p>hotomultiplier  system  approximately  10® 
away  from  the  radial  direction. 

Consequently,  the  expected  value  of  the  measured  transition  time  for  a  well-formed 
shock  was  between  0.  1  psec  and  0.  2  psec. 

Calculations  for  the  radiant  emission  from  frlly  ionized,  hydrogenic  gases 

12  13 

(Bremsstrahlung)  are  available.  '  Within  the  visible  wavelength  range  and  for  the 
anticipated  test-gas  temperatures  and  zero  magnetic  field,  the  Bremsstrahlung 
spectral  radiant  intensity  (energy/ volume/ spherical  angle/ wavelength/ time)  is  pro¬ 
portional  to  the  square  of  the  electron  density  and  is  only  a  weak  function  of  tempera¬ 
ture  (Appendix  II).  Two  wavelengths  were  chosen  from  the  high-sensitivity  range  of  the 

O  O 

photomultipliers,  namely  3800  A  (violet-ultraviolet)  and  5000  A  (green),  and  narrow 

O 

bandpass  optical  filter  systems  (half-width  <  160  A)  with  maximum  transmission  at 
each  of  these  wavelengths  were  used  in  the  experiments.  A  calibrated  incandescent 
lamp  (Fig.  8)  was  used  as  a  source  of  known  radiant  intensity,  so  that  the  absolute 
intensity  of  the  test  gas  could  be  determined  (Appendix  II). 

After  several  trial  experiments  revealed  that  the  shock  arrived  almost  simultane¬ 
ously  at  opposite  sides  of  the  shock  tube  (a  rough  check  of  azimuthal  symmetry),  the 
two  photomultipliers  were  situated  on  the  same  side  of  the  shock  tube,  at  p)ositions 
39.  5  cm  and  54.  8  cm  upstream  of  the  electrode  ring,  and  designated  PM  I  and  PM  II, 
respectively  (Fig.  8).  The  value  of  the  drive-current  maximum  usually  occurred  when 
the  shock  front  was  slightly  upstream  of  PM  I,  and  the  typical  shock  front  required 
~2  psec  to  traverse  the  distance  between  PM  1  and  PM  II.  The  variation  in  drive  current 

O 

was  ~  15  per  cent  of  the  peak  value  during  this  time.  The  3800  A  optical  filter  was 
permanently  assigned  to  PM  I,  and  the  5000  A  optical  filter  to  PM  II.  This  assignment 
was  justified  by  experiments  that  indicated  that  the  absolute  test-gas  light  signal  for 
each  optical  filter  wavelength  was  not  a  function  of  axial  position  (that  is.  of  PM  I  or 
PM  II). 

3.2  PROCEDURE 

Before  each  series  of  experiments,  the  three  oscilloscopes  (Fig.  6)  were  calibrated, 
and  each  photomultiplier  supply  voltage  was  adjusted  so  that  the  standard  incandescent 
lamp  radiation  resulted  in  an  output  signal  of  400  mv  (15  per  cent)  when  the  standard 
lamp  was  located  directly  in  front  of  the  porthole  along  the  optical  sucis  of  the 
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photomultiplier  system.  Before  each  experiment,  the  oscilloscope  scale  settings,  the 
initial  experimental  gas  pressure,  and  the  three  capacitor  bank  voltages  were  adjusted, 
and  the  values  were  recorded.  During  each  experiment,  three  oscillograms  were  taken 
(Figs.  9-17). 

The  two  photomultiplier  output  signals  were  simultaneously  displayed  on  a  dual¬ 
beam  oscilloscope,  and  the  traces  were  recorded  on  the  same  oscillogram,  the  PM  I 
trace  having  been  positioned  above  the  PM  II  trace  (Figs.  9a- 17a).  The  triggering  of 
these  traces  was  delayed  until  ~3  psec  after  the  drive -current  spark  gap  was  triggered 
(Fig.  5)  in  order  to  record  the  entire  test-gas  light  intensity  profile  for  both  PM  I  and 
PM  II.  The  sweep  speed  was  1  psec/cm  (left  to  right),  and  the  vertical  scale  was  chosen 
as  100,  200,  or  500  mv/cm  (downward  deflection  when  radiant  energy  is  received), 
varying  as  a  function  of  the  initial  experimental  gas  pressure. 

The  rate  of  change  of  drive  current  was  monitored  by  a  small-diameter  coil  of 
known  geometry  placed  between  the  copper  drive-current  plates  external  to  the  shock 
tube  (Fig.  7).  This  signal  was  integrated  (RCw  >  45)  before  being  displayed  on  a  single¬ 
beam  oscilloscope,  and  the  trace  was  recorded  on  an  oscillogram  (Figs.  9b-17b).  The 
sweep  speed  was  1  psec/cm  (left  to  right),  the  vertical  scale  was  100  mv/cm  (upward 
deflection  when  drive-current  flows),  and  the  trace  was  triggered  simultaneously  with 
the  drive-current  spark  gap  (Fig.  5). 

The  axial  magnetic  field  was  monitored  by  a  single  loop  wound  around  the  outer 
solenoid  upstream  of  PMII  (Fig.  7).  The  potential  generated  around  this  loop  was  equal 
to  the  rate  of  change  of  axial  magnetic  field  flux  through  the  loop.  Since  most  of  the 
axial  field  had  time  to  penetrate  the  steel  cylinders,  this  potential  was  approximately 
proportional  to  the  rate  of  change  of  the  axial  field  in  the  experimental  gas.  Therefore, 
this  was  designated  as  the  loop.  Data  obtained  when  examining  the  variation  of  the 
axial  field  as  a  function  of  position  indicated  that  the  maximum  axial  field  occurred  when 
the  B^  loop  signal  was  approximately  equal  to  zero.  These  data  also  indicated  that  the 
initial  value  of  the  B^  loop  signal  was  proportional  both  to  the  initial  voltage  on  the  axial 
magnetic  field  capacitor  bank  and  the  peak  value  of  the  axial  field.  Although  the  initial 
capacitor-bank  voltage  was  a  more  accurate  means  of  determining  the  peak  axial  field, 
the  B^  signal  provided  a  permanent  check  on  the  recorded  data.  The  B^  loop  signal  was 
displayed  by  the  upper  beam  of  a  dual-beam  oscilloscope  at  a  sweep  speed  of  50  psec/cm 
(left  to  right)  and  a  vertical  scale  of  either  5  or  10  v/cm,  varying  as  a  function  of  the 
initial  capacitor -bank  voltage.  The  zero  for  this  signal  was  located  on  the  middle  hori¬ 
zontal  grid.  The  oscilloscope  trace  was  triggered  simultaneously  with  the  axial-field 
spark  gap  (Fig.  5). 

The  drive-current  axial  distribution  was  monitored  by  a  loop  around  the  experi¬ 
mental  volume  in  a  plane  containing  the  shock -tube  axis  (Figs.  3  and  7).  The  potential 
generated  around  this  loop  was  equal  to  the  rate  of  change  of  azimuthal  flux,  and  the 
loop  was  designated  as  the  *  loop.  A  large  initial  *  signal  indicated  that  some  of  the 
drive  current  had  passed  from  the  inner  cylinder  to  the  outer  cylinder  at  a  point 
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upstream  of  the  electrode  ring;  consequently  the  9  signal  was  used  to  monitor  the 
initial  drive-current  breakdown  of  the  experimental  gas.  The  first  2-3  psec  of  this 
signal  were  obscured  by  hi^-frequency  noise  (>20  me),  and  a  filter  (RC  =  4  X  10  s) 
was  therefore  inserted  into  the  oscilloscope  input  line.  The  filtered  9  loop  signal  was 
displayed  by  the  lower  beam  of  a  dual-beam  oscilloscope  at  a  sweep  speed  of  1  psec/cm 
(left  to  right)  and  a  vertical  deflection  of  100  v/cm.  The  oscilloscope  trace  was 
triggered  simultaneously  with  the  drive-current  spark  gap  (Fig.  5). 

The  and  *  traces  were  recorded  on  the  same  oscillogram  (Figs.  9c-17c),  the 
former  appearing  above  the  latter. 

The  open  spark  gaps  usually  broke  down  within  3  psec  of  receiving  the  triggering 
pulse,  the  major  exception  being  the  drive-current  spark  gap  at  zero  axial  field.  This 
probably  indicates  the  importance  of  azimuthal  electron  drift  in  breaking  down  the 
experimental  gas.  In  many  experiments  run  with  zero  axial  field,  data  were  obtained 
by  using  the  •  signal  to  trigger  the  photomultiplier  oscilloscope.  No  current  trace  was 
obtained  in  these  cases  (Fig.  13). 

The  shock  tube  required  no  internal  changes  between  experiments;  and  before  each 
new  experiment  the  steady  flow  of  experimental  gas  quickly  removed  any  gaseous 
impurities  that  may  have  resulted  from  the  previous  experiment.  Each  experiment  in 
the  aeries  was  conducted  as  soon  as  the  three  capacitor  banks  had  been  recharged, 
which  required  only  a  few  moments. 

Experiments  were  run  with  the  following  mean  peak  values  of  the  axial  magnetic 
field  (in  w/m^):  0.0,  0.017,  0.034,  0.068,  0.102,  0.136,  0.170,  0.204,  0.238. 


20 


IV.  EXPERIMENTAL  RESULTS 


The  well-formed  shock  has  been  described  in  this  report  as  a  shock  for  which  the 
observed  shock-layer  total  transition  time  (based  upon  the  photomultiplier  signal  rise 
time)  is  less  than  0.2  psec.  One  of  our  outstanding  experimental  observations  is  that 
well -formed  shocks  did  not  occur  in  every  experiment.  The  experimental  results  based 
only  upon  photomultiplier  oscillogram  data  may  be  classified  as  follows: 

1 .  The  majority  (~  65  per  cent)  of  the  experiments  failed  to  produce  any  recognizable 
light  discontinuity  at  either  photomultiplier  (Fig.  9a).  These  experiments  were  almost 
always  marked  by  a  2-mc  electromagnetic  noise  that  could  not  be  associated  with  any 
known  auxiliary  circuit.  Therefore  we  concluded  that  the  2-mc  noise  was  probably 
associated  with  the  drive  current  instability  that  prevented  the  shock  from  properly 
forming. 

2.  Some  of  the  experiments  (~10  per  cent)  produced  a  well-formed  shock  at  one 
photomultiplier  and  not  at  the  other  (Figs.  10a  and  11a). 

3.  Some  of  the  experiments  (~10  per  cent)  produced  a  discernible,  well-formed 
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Fig.  9.  Data  oscillograms.  = 
0.238  w/m^,  pj  =  156pH2. 
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shock  despite  considerable  Z-mc  noise,  although  the  amount  of  the  noise  in  this  case 
was  usually  less  than  that  found  in  case  1  (Fig.  IZa). 

4.  In  the  rest  of  the  experiments  (~15  per  cent),  a  well -formed  shock  was  observed 
at  both  photomultipliers  (Figs.  13a-17a).  The  noise  level  was  extremely  low. 

The  probability  of  forming  properly  and  maintaining  a  shock  wave  in  a  given  experi¬ 
ment  appeared  to  be  a  function  of  peak  axial  magnetic  field  and  initial  experimental  gas 
pressure.  As  the  axial  magnetic  field  was  increased  from  0.  017  to  0.Z38  w/m  (and  the 
azimuthal  drift  velocity  for  electrons  diminished  by  a  factor  of  14),  the  probability  grad¬ 
ually  diminished  from  ~l/3  to  ~l/8.  For  the  case  of  zero  axial  field,  the  probability  was 
~l/6.  Furthermore,  for  any  given  value  of  axial  field,  there  appeared  to  be  a  value 
of  initial  pressure  below  which  it  was  impossible  to  produce  a  well-formed  shock.  This 
value  increased  from  ~100  p  to  ~1Z0  p  as  the  axial  field  increased  from  0.170  w/m 
to  0. 204  w/m^. 
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Fig.  13.  Data  oscillograms.  =  0.0  w/m  , 

Pj  =  llEpH^. 

These  observations  concerning  probability  indicated  the  influence  of  the  experimental 
gas  conditions  on  the  breakdown  characteristics  at  the  electrode  rings.  (Experiments 
with  zero  axial  field  established  that  Paschen's  minimum  occurred  at  a  pressure  of 
~800  ^.)  Consequently,  the  probability  of  producing  a  well-formed  shock  could  be  ten¬ 
tatively  linked  to  the  drive-current  breakdown  characteristics. 

Looked  at  as  a  whole,  the  experimental  data  represent  widely  varying  degrees  of 

success  in  obtaining  well-formed  shocks.  It  might  be  expected  that  only  the  envelope 

of  experimental  results  approached  the  theoretical  predictions,  and  it  is  upon  this 

expectation  that  the  processing  of  data  was  based.  When  at  least  20  well-formed  shocks 

were  recorded  at  the  value  of  the  axial  magnetic  field  under  investigation,  the  observed 

shock  velocity  (u  =  distance  between  successive  portholes/passage  time  of  shock  front 
s 

from  PM  I  to  PM  II)  was  compared  with  the  predicted  shock  velocity  (u^j).  The  10  shock 
waves  indicating  the  largest  percentages  of  the  theoretical  velocity  were  retained  for 
processing.  All  other  data  at  the  chosen  value  of  axial  magnetic  field  were  discarded. 
Data  obtained  at  the  maximum  axial  field  (0.238  w/m  )  constituted  the  single  exception 
to  this  procedure;  in  this  case  the  data  from  20  of  40  well -formed  shocks  were  retained 
for  processing.  The  basic  concept  was  that  those  shocks  that  represented  the  best 
shock-tube  performance  probably  adhered  most  closely  to  the  one-dimensional  theory 
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because  fewer  unexpected  influences  (for  example,  arc  instabilities)  had  acted  upon  them. 

This  procedure  required  a  method  of  computing  the  theoretical  velocity  as  a  function 
of  experimental  parameters.  When  the  axial  magnetic  field  is  zero,  one -dimensional 
theory  predicts  that  the  drive-current  sheet  will  behave  as  an  impenetrable  piston  of 

known  downstream  pressure 
be  valid.  Hence,  for  =  y. 


and  the  usual  gasdynamic  shock  relations  will 


Since 
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Fig.  16.  Data  oscillograms.  = 
0.238  w/m^,  Pj  =  136)iH2. 


Fig.  17.  Data  oscillograms.  = 
0.238  w/m^,  Pj  =  152)iH2. 


When  the  axial  magnetic  field  was  not  zero,  was  obtained  by  combining  Eq.  25 
with  the  computed  information  of  Fig.  20. 

Since  was  computed  from  measured  quantities,  the  important  errors  were  care¬ 
fully  assessed  to  allow  a  fair  comparison  of  u^j  with  as  follows: 

1.  In  Eq.  25  the  first  influence  of  the  cylindrical  geometry  upon  the  experiment  was 
found.  The  proper  choice  for  a  numerical  value  of  l/r  was  probably  between  l/r^  and 

^i  *"0 

l/r. .  The  mean  value  -= - was  used  in  the  calculations.  The  maximum  deviation 

'  1  cr.T 

10 

of  u  because  of  this  particular  choice  of  l/r,  therefore,  should  not  have  exceeded 
so 

r  -  r . 

- - — 100  =  17  per  cent. 

+  r, 
o  1 

2.  The  value  of  was  computed  from  the  measured  value  of  p^  and  the  perfect  gas 
law,  under  the  assumption  that  Tj  =  273"K.  Since  the  vacuum  gauge  was  calibrated 


25 


at  room  temperature,  the  computed  value  of  might  have  been  5  per  cent  below 

the  actual  value. 

3.  The  boundary-layer  theory  of  Appendix  I  indicates  a  14  per  cent  decrease  in  shock 

velocity  when  the  test  gas  slug  length  is  in  equilibrium. 

Since  these  errors  were  independent,  it  was  possible  to  conclude  the  fact  that  the 

expected  value  of  u„  for  zero  axial  field  should  have  been  between  0.  7  u  and  1 . 2  u 

s  so  so 

If  we  assume  that  (a)  there  was  no  axial  magnetic  field;  (b)  the  experimental  gas 
was  in  an  unexcited  state;  and  (c)  the  test  gas  was  fully  ionized,  Eqs.  21  and  22  show 


that  for  this  experiment  (hydrogen  gas,  u^ 


u, 

=  10  cm/psec),  — To.  95, 


1.5. 


Under  experimental  conditions  similar  to  those  reported  here, 


Patrick^  found  —  ~  1, 


and  Weise  et  al.^  found  — T  0. 9  and  T .  T  3  X  10^*K. 

P4  4 


Weise  et  al.  conclude  that  this 


is  the  result  of  the  delivery  of  excitation  energy  to  the  experimental  gas  by  means  other 
than  shock  heating  (for  example,  ultraviolet  radiation  from  the  drive  current)  and  in 
sufficient  amounts  to  ionize  the  experimental  gas  before  the  arrival  of  the  shock  front. 


4.  1  RAW  DATA 

a.  Photomultiplier  Oscillograms  (Figs.  I2a-17a) 

The  decrease  of  test  gas  density  (that  is,  light  intensity)  at  the  upstream  edge  of 
the  expansion  wave  was  rapid,  just  as  the  one -dimensional  theory  qualitatively  pre¬ 
dicted.  The  observed  length  of  the  test  gas  slug  was  approximately  5  cm  for  all  initial 
conditions,  which  was  consistent  with  the  boundary-layer  leakage  theory  (Appendix  I), 
but  over  far  too  wide  a  range  of  axial  fields  to  have  been  anything  more  than  a  coinci¬ 
dence.  No  precursor  light  was  observed  at  these  wavelengths.  The  success  or  failure 
of  a  given  experiment  in  producing  a  well -formed  shock  could  only  be  detected  on  the 
photomultiplier  oscillogram. 


b.  Drive-Current  Oscillograms  (Figs.  I2b-17b) 

The  drive-current  oscillograms  recorded  the  spark  gap  and  electrode  ring  break¬ 
down  delay  time  and  breakdown  noise.  The  drive  current  was  fairly  sinusoidal.  The 
variation  among  experiments  was  negligible. 

c.  Oscillograms  (Figs.  12c-17c) 

The  interrpution  of  the  signal  at  ~220  psec  was  caused  by  drive-current  break¬ 
down  noise,  and  therefore  provided  an  additional  timing  check.  Even  though  the  B^  loop 
was  in  a  plane  normal  to  the  shock-tube  axis,  and  even  though  it  was  more  than  four 
loop  diameters  distant  from  the  electrode  ring,  the  magnitude  and  duration  of  the  inter¬ 
ruption  was  proportional  to  the  2-mc  noise  found  on  the  photomultiplier  traces. 
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d.  4  Loop  (Figs.  12c-17c) 

TTie  4  loop  was  originally  designed  to  check  on  the  drive-current  breakdown,  but, 
as  shown  below,  it  yields  additional  information. 

Under  the  assumptions  that  (a)  the  drive  current  has  azimuthal  symmetry;  (b)  the 
axial  distribution  may  be  approximated  by  aJ  crossing  through  the  experimental  gas  at 
the  upstream  end  and  (l-a)J  crossing  at  the  electrode  end  (an  observation  made  after 
examining  the  disassembled  shock  tube);  and  (c)  that  the  average  drive  current  element 
moves  upstream  with  the  velocity  Uj  (pseudo-current  velocity),  Faraday's  law  reveals 
that 


aw  L  cos  wt  +  (1-0 


sin  b)t  u  cos  fa)t 


(26) 


where  is  the  magnitude  of  the  *  signal,  L  is  the  length  of  the  shock  tube,  and  w  is 
the  frequency  and  J  the  peak  value  of  the  sinusoidal  drive  current.  Thus,  at  t  =  0 


Ao  ^ 

pLJw\^r^  -r.j- 


(27) 


Since  Vj^  was  observed  to  be  <  100  v  at  t  =  0,  Eq.  27  shows  that  a  for  this  experi¬ 
ment  was  <0.03  and  may  therefore  be  neglected.  Equation  26  may  be  rearranged  with 
a  =  0  to  yield  (at  peak  current) 


J 


-  rj  j- 


(28) 


For  Vj^  :  100  V,  Uj  was  0.  70  cm/psec  for  this  experiment. 

Equation  26  plus  the  quasi-stationary  assumption  (section  2.2)  show  that  V,  is 
expected  to  be  similar  to  sin  wt  for  t  <  n/Zu,  whereas  the  oscillograms  revealed  that 
Vj^  was  proportional  to  sin  wt  for  t  <  ir/2u.  This  discrepancy  indicated  that  the  actual 
drive-current  axial  distribution  may  have  differed  from  the  prediction  of  the  one¬ 
dimensional  theory.  Although  a  4-mc  noise  was  frequently  observed  on  this  trace,  the 
2-mc  noise  was  never  observed. 


4.2  TABULATED  RESULTS 

o 

The  reduced  data  are  often  divided  into  subcritical  (0-0.  102  w/m  )  and  critical 
2 

(0.  136-0.238  w/m  )  axial  magnetic  field  groups,  in  order  to  separate  those  shocks  with¬ 
out  the  switch-on  alternative  from  those  shocks  with  the  switch-on  alternative. 

a.  Observed  Shock  Velocity  (Figs.  18-20) 

The  observed  shock  velocity  (u  )  is  compared  with  u  (not  u  , )  as  a  function  of 

s  so  S 1 

both  Pj  and  in  Figs.  18  and  19.  The  estimated  maximum  error  in  the  measurement 
of  Ug  is  10  per  cent. 
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Fig.  18.  u  and  as  functions  of  B  Fig.  19.  u  and  u_„  as  functions  of  B 

8  60  M  X  8  SO 

and  Pj,  <  0. 102  'm/vciT .  and  Pj,  B^  >  0. 136  w/m'^. 


Although  the  range  of  p^  is  too  small  to  allow  a  definite  conclusion  to  be  drawn,  the 

data  for  any  subcritical  B  indicate  that  the  dependence  of  u  upon  p ,  is  correctly  pre- 

dieted  by  Eq.  25.  The  maximum  vedues  of  u^  are  found  at  zero  B  and  are  between  80 

per  cent  and  90  per  cent  of  u^^,.  It  is  noted  that  the  predicted  range  for  u^  at  zero  B^ 

is  between  70  per  cent  and  120  per  cent  of  u  For  subcritical,  nonzero  B  ,  u  tends 

so  X  s 

to  decrease  as  B^  increases. 

The  data  for  critical  B^  exhibit  the  effect  of  the  narrowed  range  of  Pj  within  which 
well-formed  shocks  could  be  obtained.  Therefore,  no  conclusion  about  the  dependence 
of  Ug  upon  Pj  can  be  drawn  for  critical  B^. 

The  dependence  of  u^  upon  p^  is  removed  from  consideration  by  plotting  Ug/Ugj 
against  B^  (Fig.  20).  u^/u^j  does  not  approach  the  predicted  value  of  unity;  rather, 
it  rapidly  decreases  as  B^^  takes  on  a  nonzero  value  and  gradually  diminishes  as  B^^  fur¬ 
ther  increases.  The  diminishing  of  u  /u  ,  with  increasing  B  is  probably  due  to  the 

8  S 1  X 
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Fig.  20.  u  /u  ,  and  u  ./u  as  functions  of  B  . 

8  S  X  S  X  SO  X 

departure  of  actual  expansion  wave  behavior  from  the  theory.  Several  mechanisms  exist 
to  cause  the  departure,  all  of  which  are  based  upon  the  tendency  of  the  expansion-wave 
azimuthal  velocities  to  displace  the  gas  radially  and  increase  the  frequency  of  collisions 
with  the  outer  cylinder.  The  following  phenomena  result:  (a)  The  expansion  wave  gas 
properties  vary  with  radius,  constituting  a  departure  from  assumption  1  of  the  one¬ 
dimensional  theory;  (b)  the  expansion-wave  gas  is  cooled,  with  consequent  reduction 
of  conductivity,  a  departure  from  assumption  4  of  the  one-dimensional  theory;  and  (c) 
the  expansion-wave  azimuthal  velocities  will  be  lower  than  the  large  velocities  pre¬ 
dicted  by  the  one- dimensional  theory,  slowing  down  the  motion  of  the  drive  current  into 
the  test  gas.  It  might  also  be  conjectured  that  electrode  material  contaminates  the  gas 
in  the  expansion  wave,  increasing  the  mass  that  must  be  accelerated  and  decreasing 
the  resulting  shock  velocities. 

b.  Observed  Absolute  Spectral  Radiant  Intensity  (Figs.  21-23) 

The  observed  (Y^)  and  theoretical  (Y^^,  Appendix  11)  values  of  the  test-gas  absolute 
spectral  radiant  intensity  are  plotted  as  functions  of  p^  and  for  PM  I  (X  »  3800  A)  in 
Figs.  21  and  22  and  for  PM  II  (X  »  5000  A)  in  Fig.  23. 

The  PM  I  absolute  spectral  radiant  intensity  data  for  subcritical  (Fig.  21) 
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Fig.  21.  Yg  and  as  functions  of  and  Pj, 

«  0.  102  w/m^,  V  =  3800  A. 

confirm  the  assumptions  that  P^/Pj  “  4  and  that  «  3  X  10  °K  for  all  subcritical  B^. 
The  PM  I  absolute  spectral  radiant  intensity  data  for  critical  B^  (Fig.  22)  exhibit  a 
decrease  of  observed  absolute  spectral  radiant  intensity  with  increasing  B^,  but  it  is 
not  possible  to  confirm  the  switch-on  possibility  by  using  this  information  alone. 

The  PM  II  absolute  spectral  radiant-intensity  data  for  subcritical  B^  (Fig.  23)  disa¬ 
greed  with  the  theory  both  in  magnitude  and  functional  dependence  upon  p.,  Y  being 
3  2  18 
proportional  to  Pj  rather  than  Pj .  The  additional  radiation  was  possibly  due  to  line 

radiation  from  the  impurities  present  in  the  test  gas.  Consequently,  PM  II  was  used 
only  to  determine  u^. 

It  cannot  be  too  strongly  emphasized  that  the  absolute  spectral  radiant  intensity 
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Fig.  22.  Yg  and  as  functions  of  and  Pj, 
55  0.  136  w/m^.  X  =  3800  A. 


information  contained  in  the  discarded  well -formed  shock  data  was  similar  to  the 
retained  absolute  spectral  radiant-intensity  data. 

'^s/'^'se  ^sl/’^so 

The  decrease  of  absolute  spectral  radiant  intensity  at  PM  1  for  critical  B^  suggests 
the  presence  of  the  switch-on  shock  configuration.  Observed  shock-velocity  data  (Figs.  IS 
and  19)  were  combined  with  absolute  spectral  radiant-intensity  data  (Figs.  21  and 
22)  to  allow  examination  of  this  possibility.  Since  the  test-gas  absolute  spectral  radiant 
intensity  is  proportional  to  the  square  of  the  test-gas  density  (with  a  fully  ionized  test 
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Fig.  23.  Yg  and  as  functions  of  and  Pj, 
Bjj  <  0. 102  w/m^,  V  =  5000  A. 


gas  assumed),  and  the  density  ratio  is  4  for  subcritical  B^, 


Here,  Eq.  1 5  has  been  employed.  This  expression  is  especially  convenient  because 

Ysi  u 

of  the  sharp  break  in  the  functional  dependence  of -it —  upon-r —  at  u  =  2b_, . 


Fl«.  24.  Y„/Y„„  and  Y„,/Y„^  as  functions  of  b„,/u„.  \  =  3800  A. 

8  80  B 1  SO  XI  S 


Y 

Observed  values  of  y  °  and  the  theoretical  curve 
’^so 


of|ai 

so 


are  plotted  as  functions 


X 1  9 

of  — —  in  Figs.  24  and  25.  For  subcritical  B  ,  the  observed  values  of  cluster  about 

the  theoretical  curve  in  random  fashion.  The  data  for  subcritical  indicate  no  dis¬ 
cernible  functiwial  dependence  upon  u^  at  any  fixed  value  of  B^.  For  any  fixed  critical 
Bx  (Fig.  25),  however,  the  data  exhibit  a  slope  remarkably  close  to  the  theoretical  value 
of  -4.  The  scatter  at  any  fixed  critical  B^  and  b^j/u  does  not  substantially  exceed  the 
scatter  at  any  fixed  subcritical  B^  and  b^j/u  .  The  major  deviation  from  the  theory 
seems  to  be  only  the  position  of  the  knee  of  the  curve.  It  is  possible  that  the  experi¬ 
mental  points  are  as  much  as  10  per  cent  too  far  to  the  right  because:  (a)  u  is  the 
average  shock  velocity  during  the  passage  of  the  shock  from  PM  I  to  PM  II,  whereas 
the  absolute  light  intensity  was  measured  approximately  at  the  maximum  shock  velocity; 
and  (b)  any  deviation  of  the  timing  sequence  from  the  expected  finds  B  at  a  value  below 
the  maximum  value. 

In  the  particular  case  B^  b  o.  238  w/m  ,  both  a  change  of  the  experimental  pro¬ 
cedure  for  charging  the  axial -field  current  capacitor  bank,  and  the  replacement  of  a 
component  in  the  metering  section  of  the  axial-field  cap>acitor  bank  power  supply  may 
have  caused  some  additional  but  consistent  error. 

These  data  indicate  that  the  switch -on  shock  is  probably  the  stable  configuration 
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0)  Bx-O.I36w/m2  b)  B*. 0.170  w/m* 


Fig.  25.  Yg/Yg^  and  Y^j/Y^^  as  functions  of  b^j/u^  and  B^. 
?  0.  136  w/m^.  X  =  3800  A. 
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in  this  shock  tube  for  <  2b^j . 

Note  that  the  one -dimensional  shock-wave  theory  predicts*  that  the  switch-on  shock 

can  only  occur  for  >  0.453  ^  ®x  ^  although  it  was  observed 

to  occur  for  0.24.  This  discrepancy  exists  because  u^  *  u^j,  as  the  theory 

has  assumed,  but  u  ~  0.5u„, . 

8  si 

d.  Pseudo-Current  Velocity  (Fig.  26) 

The  observed  pseudo-current  velocity  (from  £q.  28)  was  found  to  diminish  with 
increasing  B  ;  therefore,  the  value  u,/u„  was  formed  and  plotted  against  B„  in  Fig.  26. 

X  J  S  X 

3  *^x4  3 

The  one -dimensional  theory  shows  that  u.,/u  ,  must  lie  between  -r-  and - +  .  The 

■'  Jl'  si  4  Ugj  4 

neglect  of  the  small  difference  between  u^j  and  u^^  results  in  the  theory  placing  the 

3  /3~  ®  1 

upper  limit  at t-  / -g-  5-— ■ ,  as  shown  in  Fig.  26.  Kemp  and  Petschek*  allow  the 

^  “yo 

precise  value  of  ^jj/u^^  to  be  computed  at  B^  =  0. 146  w/m^,  as  is  indicated  in  Fig.  26 

by  an  asterisk.  As  expected,  it  is  close  to  the  upper  limit. 

The  observed  values  of  u,/u„  are  remarkable  in  two  respects.  First,  u,/u  is  vir- 

«l  S  J  s 

tually  constant  for  all  B  .  Second,  the  average  value  of  u^u^  is  never  in  excess  of 
25  per  cent  of  •  These  two  facts  indicate  a  strong  spreading  of  the  current  sheet 

and  confirm  the  many  previous  suspicions  about  the  nature  of  the  deviation  of  the 
expansion- wave  properties  from  one -dimensional  theory. 

y 

Experiments  were  terminated  at  B^  =  0.238  w/m  because  the  absolute  spectral 
radiant  intensity  signal  was  diminishing,  and  the  range  of  useful  Pj  was  vanishing. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Magnetically  driven  shock-tube  performance  has  been  examined  from  the  standpoint 
of  two  separate  phenomena,  the  shock  layer  and  the  drive  current  carrying  expansion 
wave.  Despite  many  departures  from  the  simplified,  one -dimensional  theory,  including 
(a.)  non-steady  flow 


(b.) 


‘■o  +  ^'i 


/O 


(c.)  a  real  gas  as  the  working  fluid 

(d.)  finite  and  nonscalar  electrical  conductivity 

(e.)  viscosity  effects  in  the  boundary  layer, 

experimental  observations  indicate  that  the  one-dimensional  theory  is  adequate  for  esti¬ 
mating  the  actual  performance  in  the  annular  geometry  shock  tube. 

Shock  waves  were  found  to  exist  with  u^  <  2b^j  in  spite  of  stability  considerations 
to  the  contrary.  The  shock  transitions  were  observed,  on  the  basis  of  measured  values 
of  B^,  Pj,  Ug  and  or  p^.  to  adhere  reasonably  well  to  both  the  one -dimensional  shock- 
wave  theory  and  the  evolutionary  shock-wave  stability  theory  even  in  the  switch-on  range 
of  axial  magnetic  fields.  Well-formed  shocks  were  not  observed  in  every  experiment; 
the  cause  of  these  failures  appeared  to  be  associated  with  the  drive  current  breakdown 
characteristics. 

The  deviation  of  these  conditions  from  the  assumptions  of  one-dimensional  theory 
and  the  influence  of  the  axial  field  were  observed  from  the  relation  of  u_/u^,  and  u,/u 
to  to  influence  strongly  the  expansion  wave  behavior.  Although  Ug/Ugj  was  within 
the  expected  limits  for  zero  axial  magnetic  field,  it  continuously  decreased  as  the  axial 
field  magnitude  increased.  The  observation  that  Uj/u^  was  fairly  constant  for  all  values 
of  axial  field  but  was  less  than  one -quarter  of  the  predicted  value,  reveals  that  the  drive 
current  spreads  throughout  the  shock  tube  rather  than  concentrating  at  the  upstream 
edge  of  the  expansion  wave  as  the  one-dimensional  theory  predicts. 

The  effect  of  boundary-layer  leakage  upon  test  duration  and  shock  velocity  was  pre¬ 
dicted  on  the  basis  of  a  simple  model  and  largely  confirmed  in  these  experiments. 

Since  this  experiment  is  only  one  in  a  continuing  program  of  magnetically  driven 
shock-tube  development,  several  recommendations  are  included  to  help  direct  the  future 
effort. 

1.  If  the  operation  of  the  shock  tube  is  to  be  understood,  it  is  necessary  to  explain 
the  actual  behavior  of  the  expansion  wave.  The  complete  expansion- wave  problem  is 
extremely  complex  and  theoretically  intractable,  but  the  important  deviations  from  one¬ 
dimensional  expansion  wave  theory  might  be  explained  by  means  of  models  that  include 
only  one  additional  constraint  (for  example,  maximum  attainable  azimuthal  velocity). 

2.  In  order  to  increase  the  reliability  of  the  shock  tube  and  the  reproducibility  of 
the  experimental  observations,  arc  theory  should  be  considered  in  an  effort  to  improve 
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the  breakdown  behavior.  It  is  possible  that  the  substitution  of  a  different  experimental 
gas  would  be  an  improvement,  even  though  the  shock  velocity  decreased. 

3.  The  observed  test  duration  was  quite  small  (~0.  5  psec),  and  few  experiments 
can  be  completed  within  this  period.  Theoretically,  the  maximum  test  duration  for  this 
experiment  is  ~2  psec,  and  attempts  should  be  made  to  approach  this  value. 

4.  The  knowledge  of  any  local  test  gas  and  expansion- wave  properties  would  con¬ 
tribute  significantly  to  an  understanding  of  the  shock-tube  performance.  Devices  that 
indicate  any  possibility  of  success  in  this  application  should  be  thoroughly  investigated. 

5.  Complicated  flow  patterns  can  often  be  deciphered  from  smear  photographs.  The 
current  photographic  technology  is  such  that  smear  photographs  can  be  made  of  the  test 
gas  in  this  experiment,  and  the  possibility  of  obtaining  such  information  should  not  be 
overlooked  in  future  experiments. 
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APPENDIX  1 


ESTIMATES  OF  THE  TEST  DURATION 
1.  MINIMUM  TEST  DURATION 

2 

The  drive-current  pressure  is  proportional  to  l/r  (Eq.  24),  causing  the  upstream 
edge  of  the  expansion  wave  to  bend  downstream  with  increasing  radius,  with  the  result 
that  a  finite  volume  of  test  gas  exists  even  when  the  drive  current  overtakes  part  of  the 


using  the  nomenclature  of  Fig.  27  and  assuming^; 


(i)  steady  flow  (shock  front  and  expansion  wave  move  with  equal  velocities); 

(ii)  an  impenetrable  drive-current  interface  (the  flow  direction  is  parallel  to  the 
intei'face); 


(iii)  L(r)  is  small  enough  to  insure  that  all  properties  are  functions  of  r  alone; 

(iv)  the  axial  magnetic  field  is  zero; 

(v)  p  . 


k,  +  1 

k,  -  1 


r 


(vi)  Piu^=^: 

(vii)  =  L 

(viii)  may  have  any  value  at  r 

as  a  sink  for  the  gas); 


-  —  -  0  at  r  =  r^  (the  drive  current  has  overtaken  the  shock  front); 

r^  (Appendix  1.  2  shows  that  the  outer  wall  behaves 


(ix)  — 1. 


Defining 


and 

P  "  P4''4l'rU4  -  cu_^,  (31) 


we  have  the  conservation  equations 


d  I  \ 


Pl^ir 


(continuity) 


(32) 
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^  (P)  =  PjUjr  J  (x-momentum) 

^  (av^)  =  Piu5|“  lJ  ►  (r-momentum) 

Equations  32  and  33  can  be  integrated  directly  to  yield 

Since  =  p/a  (Eq.  31),  and  r/r^  —  1,  it  follows  that 


(33) 


(34) 


(35) 


(36) 


“4  =“,  «( 


2  In 


(^) 


-1 


(37) 


and  also  that 
P  »  P, 


(38) 


By  neglecting  ^ 
be  shown  to  be  valid  for  r/r^  >  1.1),  (34)  becomes 

k2  , / _  \ 2 


*’i“4 

comparison  with-— — in  Eq.  34  (an  approximation  which  can 
rv4 


which  may  be  integrated  to  yield 


3/2 


av^  «. 


whence 


Also 


2  2  /  r 


1/2 


(39) 


(40) 


(41) 


(42) 


Note  that  the  gas  properties  enter  into  the  equation  only  for  L(r).  Since  dxj/dr 

U4/V4. 
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Equation  43  may  be  integrated  to  yield 
.3/2 


(44) 


This  equation  is  compared  below  with  the  data  of  Keck,  Fishman,  and  Petschek^ 
as  a  rough  check  of  the  validity  of  the  assumptions  used  in  the  derivations  above. 

r 

X 1 

-^(Eq.  44) 

X.  (. 

-p—  (Keck,  et  al.  ) 

1 

0 

0 

1.5 

0.29 

0.2  -  0.4 

2.0 

0.  82 

0.  6  -  0. 9 

2.5 

1.  50 

1.2  -  2.0 

The  agreement  between  the  ccKnputed  and  observed  values  is  excellent,  whence  it 
may  be  concluded  that  this  model  properly  describes  the  flow  pattern  when  the  drive 
current  overtakes  the  shock  front,  because  of  the  variation  of  pressure  with  radius.  The 
other  important  properties  are  presented  below  for  k^  =  5/3  and  for  the  range  of  values 
of  r/r.  found  in  the  shock  tube  investigated  here.  They  are  valid  only  for  r/r.  >  1.1, 

^  /  r.  \  *■.  u . 

because  of  the  assumption  that  "JT 

^  (Eq.  44) 
*^i  *^i 

^(Eq.  42) 

X  X 

-  ^(shock  front) 

1  0 

0 

0 

1.1  0.03 

0.  10 

-0.  07 

1.2  0.07 

0.  14 

-0.07 

1.3  0.13 

0.  17 

-0.04 

1.4  0.20 

0.20 

-0.00 

The  predicted  values  of  x^/ r^  in  the  range  of  validity  of  r/ r^  indicate  that  the  shock 
front  is  tilted  slightly  downstream  with  increasing  radius.  The  average  angle  of  shock- 
front  tilt  between  r/r^=  1.1  and  r/r^  =  1.4  is  approximately  10”  (a  value  that  is 

independent  of  u^).  The  time  required  for  the  shock  front  to  pass  a  given  axial  point 
would  be  approximately  0.  07  r^u^  <  0.  1  jis  for  this  experiment.  The  minimum  length 
of  the  test-gas  volume  at  the  mean  radius  (r/r.  =  1. 2)  is  taken  to  be  the  appropriate 
measure  of  the  minimum  length  of  the  test-gas  slug,  and  this  length  has  the  value 
“0.7  cm  for  this  experiment. 


2.  MAXIMUM  TEST  DURATION 
8  9 

Roshko  and  Duff  have  shown  that  the  transfer  of  test  gas  into  the  boundary  layer 
prevents  the  test-gas  volume  in  combustion-driven  shock  tubes  from  increasing  indefi¬ 
nitely.  Since  the  test  gas  that  condenses  on  the  wall  of  the  magnetically  driven  shock 
tube  suffers  a  large  decrease  in  conductivity,  the  condensed  gas  therefore  leaks  past 
the  drive  current  without  difficulty,  and  a  phenomenon  similar  to  that  observed  in 
combustion-driven  shock  tubes  can  result.  Following  Roshko,  (using  the  nomenclature 
of  Fig.  27).  but  limiting  the  derivation  to  strong  gasdynamic  shocks  (M^  =  w)  in  hydro¬ 
gen  and  neglecting  any  boundary  layer  mass  transfer  at  r  =  r .  (assuming  the 
drive-current  configuration  of  Appendix  1. 1).  we  can  make  an  estimate  of  the  maximum 
test-gas  volume. 

Since 


=  2’'*'o  yp4W^w“l^o*  (^5) 

where  m^,  is  the  rate  of  mass  transfer  out  of  the  test-gas  volume,  p  is  the  dynamic  vis¬ 
cosity,  and  W  refers  to  conditions  at  the  outer  wall  of  the  shock  tube,  and 


/  ^  4W 

^‘4W  ^1  y  Tj  ■■  ^1 

(approximate  Sutherland  law^^) 

(46) 

P4'r4  '‘Pl'^4 

»’4W  -  T,^  T, 

(perfect  gas) 

(47) 

2 

T  « 4^ 

4  16  R^ 

(Eq.  13). 

(48) 

it  follows  therefore  that 


O 


m 


c 


y  «4Ti 


(49) 


The  rate  of  mass  transfer  into  the  test-gas  volume  is  referred  to  as  m^,  and  is  given 
by  the  relation 


(50) 
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At  equilibrium,  dL^/dt  =  0,  =  ^oe'  "  '”c’ 


L 

oe 


(51) 


For  this  experiment.  «  4  cm.  Since  the  value  of  »  4  cm  is  somewhat  larger 
than  the  value  obtained  in  section  1  for  the  minimum  length  of  the  test-gas  volume 
(L^  «  0.  7  cm),  the  conclusion  is  therefore  drawn  that  the  shock  wave  separates  from 
the  drive-current  interface,  but  that  the  shock  phenomenon  retains  some  of,the  features 
predicted  above  (e.  g. ,  the  downstream  tilting  of  the  shock  front  and  current  interface 
with  increasing  radius). 

O  o 

The  amount  of  time  required  for  equilibrium  between  m^  and  m^  to  be  established 
(under  the  assumption  of  a  step-function  change  in  the  shock  velocity)  is  given  by 


6L 


O®  ~  o  c 

- *2.5  (is. 


(52) 


where  >  0.  7  for  t  >  t^.  Consequently,  equilibrium  should  be  established  by  the 
time  the  shock  front  reaches  PMI  (t~5  ps),  and  the  anticipated  experimental  value  of 
the  test-gas  slug  length  is  therefore  ~4  cm. 

The  removal  of  test  gas  by  the  boundary  layer  causes  the  drive-current  interface 
to  accelerate  and  the  shock  to  decelerate  until  they  travel  at  equal  velocities.  The  equi¬ 
librium  value  of  this  velocity  is  estimated  by  means  of  the  simple  model  described 
below. 

In  a  reference  frame  attached  to  the  shock  front  and  current  interface,  the  conser¬ 
vation  of  axial  momentum  at  equilibrium  is  expressed  by  the  relation 


2  /  2  2  \ 

jUj  -  r^  j  +  skin  friction  that  is  due  to  moving  wall 

T(r^  -  r?  1+  x-momentum  outflow  associated  with  m 
\o  I  )  c 


Pa  > 


(53) 


In  this  reference  frame  the  wall  velocity  is  equal  to  u^,  and  is  greater  than  4u^^.  Con¬ 
sequently,  the  test  gas  appears  to  be  standing  still,  and  any  x-momentum  imparted  to 
the  test  gas  in  the  boundary  layer  appears  to  be  the  result  of  the  skin  friction,  that  is. 


skin  friction  *  outflow  associated  with  m  . 

c 

Therefore 


Since  the  ideal  gasdynamic  relation  (Eq.  23)  is 


(54) 


(55) 
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it  may  be  concluded  that  the  loss  of  test  gas  to  the  boundary  layer  results  in  a  14  per 

cent  decrease  of  the  observed  shock  velocity  from  the  theoretical  shock  velocity  based 

upon  p.  and  p,.  This  14  per  cent  decrease  in  the  observed  shock  velocity  is  apparent 
^  1  () 
in  the  data  of  experiments  conducted  by  Duff. 
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APPENDIX  II 


ABSOLUTE  SPECTRAL  RADIANT-INTENSITY  INFORMATION 


I.  AB;jOLUTE  spectral  radiant  intensity  of  fully  IONIZED 
HYDROGEN  FOR  ZERO  MAGNETIC  FIELD 
12 

Janes  and  Koritz  numerically  calculated  the  absolute  Bremsstrahlung  spectral 
radiant  intensity  (Y^^  =  energy/volume/spherical  angle/wavelength/time)  for  the  tem¬ 
perature  and  wavelength  ranges  of  this  experiment  and  zero  magnetic  field,  and  the  cal¬ 
culations  are  reproduced  on  Fig.  28  (X.  =  3800  A)  and  Fig.  29  (X  =  5000  A).  The  computed 
results  clearly  indicate  that  Y^^  is  proportional  to  the  test-gas  electron  number  density 
(n)  squared,  and  that  Y^^  is  only  a  weak  function  of  the  test-gas  temperature  (T).  In 

fact,  assuming  that  T  =  3  X  10^*K,^  we  find  that  the  value  Y  /n^  =  3.1  X  10  for 
o  4  5 

X  =  3800  A  is  accurate  within  10  per  cent  for  5x10  <T<5XlO,  and  the  value 

Y  /n^  =  2.0  X  10  for  X  =  5000  A  is  accurate  within  10  per  cent  for  5X  10^<  T  <4X  10^. 

Oster  has  derived  the  following  analytical  expressions  for  Y  /n  ,  based  upon 

SO 

the  assumption  that  the  radiating  electrons  move  in  straight  lines; 


*so 

a 

In  (pXT^/^) 

(classical  range) 

(57) 

2 

n 

'  x^nTt 

^so 

a 

In  (yXT) 

(quantum  range) 

(58) 

n2 

Here,  a  =  0.90  X  10  p  =  1.65  X  10  and  y  =  1.55  X  10^^.  The  classical  range  is 
2  2 

T  ^  (y/P)  »  and  the  quantum  range  is  T  ^(y/P)  •  Both  expressions  are  valid  only  for 
T  »  1.5/(X)^/^  ~  3  X  10^'K. 

Referring  to  Figs.  28  and  29,  it  is  apparent  that  the  computed  and  analytical  methods 
2  5 

of  computing  Y^^/n  agree  for  T  >  3  X  10  ’K.  The  analytical  formulation  has  the  addi¬ 
tional  advantage  of  allowing  T  and  n^  to  be  directly  computed  from  the  values  of  Y^  at 
two  separate  wavelengths  (I  and  II),  as  follows: 


(classical  range) 


(quantum  range) 


(classical  range) 


(59) 


(60) 


(61) 
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Fig.  28.  as  a  function  of  T,  \  =  3800  A. 


Fig.  29.  Yg^/n^  as  a  function  of  T,  X  =  5000  A. 
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(quantum  range) 


(62) 


where 


r  = 


Y 

V 


2.  COMPUTATION  OF  THE  OBSERVED  ABSOLUTE  SPECTRAL 
RADIANT  INTENSITY 


The  absolute  spectral  radiant  intensity  of  the  test  gas  was  calculated  by  comparing 
the  observed  photomultiplier  signal  with  the  signal  that  is  due  to  a  source  of  known  inten¬ 
sity,  as  described  below. 

14 

Barbrow  has  described  a  method  for  computing,  with  5  per  cent  accuracy,  the 
spectral  radiant  intensity  (W  =  energy /spherical  angle/wavelength/time)  of  a  tungsten 

filament  incandescent  lamp  of  known  radiant  flux  (P  =  energy/time)  and  approximate 

15 

filament  temperature.  Forsythe  provides  the  data  that  establish  the  facts  that 


W  ,  =  3.6  X  10^  P 
cl 

W  „  =  3.6  X  10^  P, 
cll 

o  o 

where  c  refers  to  the  calibration  value,  Xj  =  3800  A,  and  Xjj  =  5000  A. 
Neglecting  reflection,  refraction,  etc.,  we  have 


W  =  Y  AV, 

so  so 


(63) 


(64) 


where  AV  is  a  small  volume  of  test  gas  distant  from  and  visible  to  the  photomultiplier. 

Since  W^.An^,  where  V  is  the  photomultiplier  output  signal,  and 

Ail  is  the  spherical  angle  of  radiation  intercepted  by  the  photomultiplier  optical  system, 
it  follows  that 


Y  AVa« 

so  so 

W  An 

c  c 


(65) 


Under  the  assumptions  of  a  compression  ratio  of  4  across  the  shock  and  complete  ioni¬ 
zation  of  the  test  gas,  n^  =  8.0  X  10^°  Pj  (pj  in  p),^^  and  therefore  from  Appendix  II  it 
follows  that 
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(66) 


V.  ■  x  A 

The  method  of  calibration  (section  3.2)  fixes  V  b  0.40  and  ^(2  =  AO  Since  P 

-R  -a  ^  c  so 

AVj  =  0.91  X  10  ,  and  AVjj  =  1.54  X  10  ,  it  follows  that 

''sol  =  X  ‘O'*  A 

(Pj  in  ti) 

''.oU-'‘’x  10-‘p5 

If  the  observed  photomultiplier  signal  magnitude  is  V^,  then 


Except  when  otherwise  noted,  the  units  above  are  mks. 


=  0.94. 


(67) 


(68) 
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APPENDIX  III 


AUXIUARY  EQUIPMENT 

Vacuum  Equipment 

1  National  Research  Corporation  H-2-SP  diffusion  pump 
1  National  Research  Corporation  MO-76390  molecular  trap 
1  National  Research  Corporation  520  ionization  gauge 
1  W.  M.  Welch  Manufacturing  Company  1402B  mechanical  pump 

1  Vacuum  Engineering  Corporation  VL  variable  leak  valve 

Power  Supplies 

2  New  Jersey  Engineering  Corporation  H-10-50  power  supplies 
1  Potter  Company  PHV30-IM60V  power  supply 

Capacitance 

43  Aerovox  Corporation  PX50J5  capacitors 

12  Aerovox  Corporation  PX50D4  capacitors 

4  Cornell -Dubilier  Electronics  30-DK-T5  capacitors 

Triggering  Equipment 

1  Dumont  Corporation  404  pulse  generator 

3  Avco  Corporation  designed  1.  5  KV  pulse  generators 

1  M.I.T.  designed  100*300  psec  delay  circuit 

3  Laboratory  for  Electronics,  Inc.,  T-71 358-2  pulse  transformers 
Photomultiplier  Equipment 

2  Radio  Corporation  of  America  1P21  photomultipliers 

2  Northeast  Scientific  Corporation  RE- 1602  power  supplies 
2  Avco  designed  photomultiplier  circuits 

1  Bausch  and  Lomb,  Inc.,  33-78-38  second-order  interference  filter 

1  Bausch  and  Lomb,  Inc.,  33-78-50  second-order  interference  filter 

2  Corning  Glass  Works  4308  cutoff  filters 

1  Electrical  Testing  Laboratories,  Inc.,  8103  standard  lamp 

Oscilloscopes 


1 

Tektronix,  Inc., 

555  dual-beam  oscilloscope 

2 

Tektronix,  Inc., 

Type  G  preamplifiers 

1 

Tektronix,  Inc., 

551  dual-beam  oscilloscope 

2 

Tektronix,  Inc., 

Type  L  preamplifiers 

1 

Tektronix,  Inc., 

535  single-beam  oscilloscope 

1 

Tektronix,  Inc., 

Type  B  preamplifier 
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